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C§ . ABSTRACT 



^ ■ We undertook a long term project, DIRECT, to obtain the direct distances to two 



important galaxies in the cosmological distance ladder - M31 and M33, using detached 
eclipsing binaries (DEBs) and Cepheids. While rare and difficult to detect, detached 
eclipsing binaries provide us with the potential to determine these distances with an 
accuracy better than 5%. The massive photometry obtained in order to detect DEBs 
provides us with the light curves for the Cepheid variables. These are essential to the 
parallel project to derive direct Baade-Wesselink distances to Cepheids in M31 and 
M33. For both Cepheids and eclipsing binaries the distance estimates will be free of 
any intermediate steps. 

As a first step of the DIRECT project, between September 1996 and January 1997 
we have obtained 36 full nights on the Michigan-Dartmouth-MIT (MDM) 1.3-meter 



^Based on the observations collected at the Michigan-Dartmouth-MIT (MDM) 1.3-meter telescope and at the 
F. L. Whipple Observatory (FLWO) 1.2-meter telescope 

^On leave from N. Copernicus Astronomical Center, Bartycka 18, Warszawa 00-716, Poland 
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telescope and 45 full/partial nights on the F. L. Whipple Observatory (FLWO) 
1.2-meter telescope to search for detached eclipsing binaries and new Cepheids in 
the M31 and the M33 galaxies. In this paper, second in the series, we present the 
catalog of variable stars, most of them newly detected, found in the field M31A 
(02000.0, <J2000 = 11.34 deg, 41.73 deg). We have found 75 variable stars: 15 eclipsing 
binaries, 43 Cepheids and 17 other periodic, possible long period or non-periodic 
variables. The catalog of variables, as well as their photometry and finding charts, are 
available using the anonjrmous ftp service and the WWW. The CCD frames are available 
on request. 

Subject headings: distance scale — galaxies:individual(M31,M33) — eclipsing binaries — 
Cepheids 

1. Introduction 

The two nearby galaxies - M31 and M33, are stepping stones to most of our current effort 
to understand the evolving universe at large scales. First, they are essential to the calibration of 
the extragalactic distance scale (Jacoby et al. 1992; Tonry et al. 1997). Second, they constrain 
population synthesis models for early galaxy formation and evolution, and provide the stellar 
luminosity calibration. There is one simple requirement for all this - accurate distances. 

Detached eclipsing binaries (DEBs) have the potential to establish distances to M31 and 
M33 with an unprecedented accuracy of better than 5% and possibly to better than 1%. These 
distances are now known to no better than 10-15%, as there are discrepancies of 0.2 — 0.3 mag 
between RR Lyrae and Cepheids distance indicators (e.g. Huterer, Sasselov & Schcchter 1995). 
Detached eclipsing binaries (for reviews see Andersen 1991, Paczyhski 1997) offer a single step 
distance determination to nearby galaxies and may therefore provide an accurate zero point 
calibration - a major step towards very accurate determination of the Hubble constant, presently 
an important but daunting problem for astrophysicists (see the papers from the recent "Debate 
on the Scale of the Universe": Tammann 1996, van den Bergh 1996). 

The detached eclipsing binaries have yet to be used (Huterer et al. 1995; Hilditch 1996) as 
distance indicators to M31 and M33. According to Hilditch (1996), there are about 60 eclipsing 
binaries of all kinds known in M31 (Gaposchkin 1962; Baade &: Swope 1963; Baade &: Swope 1965) 
and only one in M33 (Hubble 1929)! Only now does the availability of large format CCD detectors 
and inexpensive CPUs make it possible to organize a massive search for periodic variables, which 
will produce a handful of good DEB candidates. These can then be spectroscopically foUowed-up 
with the powerful new 6.5-10 meter telescopes. 

The study of Cepheids in M31 and M33 has a venerable history (Hubble 1926, 1929; 
Gaposchkin 1962; Baade &; Swope 1963; Baade & Swope 1965). In the 80's Freedman &; 
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Madore (1990) and Preedman, Wilson &; Madore (1991) studied small samples of the earlier 

discovered Cepheids, to build PL relations in M31 and M33, respectively. However, both the 
sparse photometry and the small samples do not provide a good basis for obtaining direct 
Baade-Wesselink distances (e.g. Krockenberger, Sasselov & Noyes 1997) to Cepheids - the need 
for new digital photometry has been long overdue. Recently, Magnier et al. (1997) surveyed large 
portions of M31, which have previously been ignored, and found some 130 new Cepheid variable 
candidates. Their light curves are however rather sparsely sampled and in V band only. 

In Kaluzny et al. (1998) (hereafter: Paper I), the first paper of the series, we presented the 
catalog of variable stars found in one of the fields in M31, called M31B. Here we present the 
catalog of variables from the neighboring field M31A. In Scc.2 wc discuss the selection of the fields 
in M31 and the observations. In Sec. 3 we describe the data reduction and calibration. In Sec. 4 we 
discuss briefly the automatic selection we used for finding the variable stars. In Sec. 5 we discuss 
the classification of the variables. In Sec.6 we present the catalog of variable stars. 

2. Fields selection and observations 

M31 was primarily observed with the McGraw-Hill 1.3-meter telescope at the MDM 
Observatory. We used the front-illuminated, Loral 2048^ CCD Wilbur (Metzger, Tonry & Luppino 
1993), which at the //7.5 station of the 1.3-meter has a pixel scale of 0.32 arcsec/pixel and field 
of view of roughly 11 arcmin. We used Kitt Peak Johnson-Cousins BVI filters. Some data for 
M31 were also obtained with the 1.2-meter telescope at the FLWO, where we used "AndyCam" 
with thinned, back-side illuminated, AR coated Loral 2048^ CCD. The pixel scale happens to be 
essentially the same as at the MDM 1.3-meter telescope. We used standard Johnson- Cousins BVI 
filters. 

Fields in M31 were selected using the MIT photometric survey of M31 by Magnier et 
al. (1992) and Haiman et al. (1994) (see Paper I, Fig.l). We selected six 11' x 11' fields, M31A-F, 
four of them (A-D) concentrated on the rich spiral arm in the north-eastern part of M31, one (E) 
coinciding with the region of M31 searched for microlensing by Crotts & Tomaney (1996), and one 
(F) containing the giant star formation region known as NGC206 (observed by Baade & Swope 
1963). Fields A-C were observed during September and October 1996 5-8 times per night in the 
V band, resulting in total of 110-160 V exposures per field. Fields D-F were observed once a night 
in the F-band. Some exposures in B and / bands were also taken. M31 was also occasionally 
observed at the FLWO 1.2- meter telescope, whose main target was M33. 

In this paper we present the results for the M31A field. We obtained for this field useful data 
during 29 nights at the MDM, collecting a total of 109 exposures of 900 sec in V, 27 exposures of 
600 sec in I and 2 exposures of 1200 sec in B. We also obtained for this field useful data during 15 
nights at the FLWO, collecting a total of 8 exposures of 900 sec in V and 18 exposures of 600 sec 
in /. The complete list of exposures for this field and related data files are available through 
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anonymous ftp on cfa-ftp.harvard.edu, in pub/kstanek/DIRECT directory. Please retrieve the 
README file for instructions. Additional information on the DIRECT project is available through 



the WWW at tittpV/cf a-www.harvard.edu/~kstanek/DIRECT/ 



3. Data reduction, calibration and astrometry 

The details of the reduction procedure are given in Paper I. Preliminary processing of the 
CCD frames was done with the standard routines in the IRAF-CCDPROC package.^ Stellar profile 
photometry was extracted using the Daophot/Allstar package (Stetson 1987, 1991). We selected a 
"template" frame for each filter using a single frame of particularly good quality. These template 
images were reduced in a standard way (Paper I). Other images were reduced using Allstar in the 
fixed-position-mode using as an input the transformed object list from the template frames. For 
each frame the list of instrumental photometry derived for a given frame was transformed to the 
common instrumental system of the appropriate "template" image. Photometry obtained for the 

V I filters was combined into separate data bases. Unlike for the M31B field, M31A images 
obtained at the FLWO were reduced using MDM "templates". 

To obtain the photometric calibration we observed 4 Landolt (1992) fields containing a total 
of 18 standards stars. These fields were observed through BVI filters at air-masses ranging from 
1.2 to 1.70. The transformation from the instrumental to the standard system was derived and 
described in Paper I. The derived transformation satisfactorily reproduces the V magnitudes and 

V — I colors. The B — V transformation reproduces the standard system poorly, and we decided 
to drop the B data from our analysis, especially since we took only 2 B frames. 

To check the internal consistency of our photometry we compared the photometry for 20 
y < 20 and 47 / < 20 common stars in the overlap region between the fields M31A and M31B 
(Paper I, Fig.l). There was an offset of 0.022 mag in V and 0.018 mag in /, i.e. well within 
our estimate of the 0.05 mag systematic error discussed in Paper I. We also derived equatorial 
coordinates for all objects included in the data bases for the V filter. The transformation from 
rectangular coordinates to equatorial coordinates was derived using ~ 200 stars identified in the 
list published by Magnier et al. (1992). 



4. Selection of variables 

The procedure for selecting the variables was described in detail in Paper I, so here we only 
give a short description, noting changes when necessary. The reduction procedure described in 
previous section produces databases of calibrated V and I magnitudes and their standard errors. 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Associations 
of Universities for Research in Astronomy, Inc., under cooperative agreement with the NSF 



- 5 - 

T , , I , I I I , , I I I I r 




Fig. 1. — Distributions in V (dashed line) and / (continuous line) of stars in the field M31A. 

The V database for M31A field contains 10084 stars with up to 117 measurements, and the / 
database contains 21341 stars with up to 45 measurements. Fig.|^ shows the distributions of stars 
as a function of mean V oi I magnitude. As can be seen from the shape of the histograms, our 
completeness starts to drop rapidly at about V ^ 22 and I ~ 20.5. The primary reason for this 
difference in the depth of the photometry between V and / is the level of the combined sky and 
background light, which is about three times higher in the I filter than in the V filter. 

The measurements flagged as "bad" and measurements with errors exceeding the average 
error by more than 4cr are removed (Paper I). Usually 0~10 points are removed, leaving the 
majority of stars with roughly Ngood ~ 105 — 117 V measurements. For further analysis we use 
only those stars which have at least Ngood > Nmax/2 (= 58) measurements. There are 8521 such 
stars in the V database of the M31A field. 

Our next goal is to select objectively a sample of variable stars from the total sample defined 
above. There are many ways to proceed, and we largely follow the approach of Stetson (1996). 
The procedure is described in more detail in Paper I. In short, for each star we compute the 
Stetson's variability index Js (Paper I, Eq.7), and stars with values exceeding some minimum 
value Js,min are considered candidate variables. The definition of Stetson's variability index 
includes the standard errors of individual observations. If, for some reason, these errors were over- 
or underestimated, we would either miss real variables, or select spurious variables as real ones. 
Using the procedure described in Paper I, we scale the Daophot errors to better represent the 
"true" photometric errors. We then select the candidate variable stars by computing the value of 
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Fig. 2. — Variability index Js vs. mean V magnitude for 8521 stars in the field M31A with 
Ngood > 58. Dashed line at Js = 0.75 defines the cutoff applied for variability. 



Js for the stars in our V database. We used a cutoff of Js,min = 0.75 and additional cuts described 
in Paper I to select 183 candidate variable stars (about 2% of the total number of 8521). In Fig.^ 
we plot the variability index Js vs. apparent visual magnitude V for 8521 stars with Ngood > 58. 



5. Period determination, classification of variables 

We based our candidate variables selection on the V band data collected at the MDM and 
the FLWO telescopes. We also have the I band data for the field, up to 27 MDM epochs and up 
to 18 FLWO epochs. As mentioned above, the / photometry is not as deep as the V photometry, 
so some of the candidate variable stars do not have an / counterpart. We will therefore not use 
the / data for the period determination and broad classification of the variables. We will however 
use the / data for the "final" classification of some variables. 

Next we searched for the periodicities for all 183 candidate variables, using a variant of the 
Lafler-Kinman (1965) technique proposed by Stetson (1996). Starting with the minimum period 
of 0.25 days, successive trial periods are chosen so 

P-I -P-^-^ (I) 



where At = t]\j — ti is the time span of the series. The maximum period considered is 
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At = 55.8 days. For each candidate variable 10 best trial periods are selected (Paper I) and then 
used in our classification scheme. 

The variables we are most interested in are Cepheids and eclipsing binaries (EBs). We 

therefore searched our sample of variable stars for these two classes of variables. As mentioned 
before, for the broad classification of variables wc restricted ourselves to the V band data. We will, 
however, present and use the / band data, when available, when discussing some of the individual 
variable stars. 

In the search for Cepheids we followed the approach by Stetson (1996) of fitting template light 
curves to the data. We used the parameterization of Cepheid light curves in the V band as given 
by Stetson (1996). Unlike for the M31B field, we classified the star as a Cepheid if the reduced 
X^/^DOF of the fit was factor of 3 smaller than the reduced x^/^dof of a straight line fit, not a 
factor of 2 smaller. We also allowed for periods longer than 3 days, not longer than 4 days as in 
Paper I, which resulted in finding two possible Cepheids with periods 3 < P < 4 days. There was 
a total of 51 variables passing all of the criteria. Their parameters and light curves are presented 
in the Sections 6.2, 6.3. 

For eclipsing binaries (EBs) we used similar search strategy to that described in detail in 
Paper I, but we simplified the condition for the reduced x^/^dof of the fit to be at least 1.75 
smaller than the reduced X? /^dof of a straight line fit. Within our assumption the light curve of 
an EB is determined by nine parameters: the period, the zero point of the phase, the eccentricity, 
the longitude of periastron, the radii of the two stars relative to the binary separation, the 
inclination angle, the fraction of light coming from the bigger star and the uneclipsed magnitude. 
A total of 15 variables passing all of the criteria and their parameters and light curves are 
presented in the Section 6.1. 

After wc selected 15 eclipsing binaries and 51 possible Cepheids, we were left with 118 "other" 
variable stars. After raising the threshold of the variability index to Js,min = 1-2 (Paper I) we are 
left with 29 variables, which we preliminary classify as "miscellaneous" . We then go back to the 
CCD frames and try to see by eye if the inferred variability is indeed there, especially in cases 
when the light curve is very noisy/chaotic. We decided to remove 2 dubious Cepheids and 18 
dubious miscellaneous variables from the sample, which leaves 11 variables which we classify as 
miscellaneous. Their parameters and light curves are presented in the Section 6.4. 



6. Catalog of variables 

In this section wc present light curves and some discussion of the 75 variable stars 
discovered in our survey. Complete V and (when available) / photometry and 128 x 128 pixel 
(~ 40" X 40") V finding charts for all variables are available through the anonjrmous ftp on 
cfa-ftp.harvard.edu, in pub/kstanek/DIRECT directory. Please retrieve the README file for 
the instructions and the list of files. These data can also be accessed through the WWW at the 
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littp://cf a-www.harvard.edu/~kstanek/DIRECT/. The CCD frames for this field, and also for 



the M31B field, are available on request from K. Z. Stanek (kstanek@cfa.harvard.edu). 

The variable stars are named according to the following convention: letter V for "variable", 
the number of the star in the V database, then the letter "D" for our project, DIRECT, followed 
by the name of the field, in this case (M)31A, e.g. V3407 D31A. Tables |l], |, | and | list the 
variable stars sorted broadly by four categories: eclipsing binaries, Cepheids, other periodic 
variables and "miscellaneous" variables, in our case meaning "variables with no clear periodicity". 
Some of the variables which were found independently by survey of Magnier et al. (1997) are 
denoted in the "Comments" by "Ma97 ID" , where the "ID" is the identification number assigned 
by Magnier at al. (1997). We also identify several variables found by us in Paper I. 

Note that this is a first step in a long-term project and we are planning to collect additional 
data and information of various kind for this and other fields. As a result, the current catalog 
might undergo changes, due to addition, deletion or re-classification of some variables. 



6.1. Eclipsing binaries 

In Table [l| we present the parameters of the 15 eclipsing binaries in the M31A field. The 
lightcurves of these variables are shown in Figs.^-^, along with the simple eclipsing binary models 
discussed in the Paper I. The variables are sorted in the Table |l] by the increasing value of the 
period P. For each eclipsing binary we present its name, 2000.0 coordinates (in degrees), value of 
the variability index Js, period P, magnitudes Vmax and Imax of the system outside of the eclipse, 
and the radii of the binary components i?2 in the units of the orbital separation. We also 
give the inclination angle of the binary orbit to the line of sight i and the excentricity of the orbit 
e. The reader should bear in mind that the values of Vmax, Imax, Ri, R2, i and e are derived 
with a straightforward model of the eclipsing system, so they should be treated only as reasonable 
estimates of the "true" value. 

One of the eclipsing binaries found, V4636 D31A, is a good DEB candidate, with reasonably 
deep eclipses and the ellipticity indicating that the system is young and unevolved. However, 
much better light curve is necessary to accurately establish the properties of the system. 



6.2. Cepheids 

In Table || we present the parameters of 43 Cepheids in the M31A field, sorted by the period 
P. For each Cepheid we present its name, 2000.0 coordinates, value of the variability index 
Js, period P, flux-weighted mean magnitudes {V) and (when available) (/), and the F-band 
amplitude of the variation A. In Figs.^-^ we show the phased V, I lightcurves of our Cepheids. 
Also shown is the best fit template lightcurve (Stetson 1996), which was fitted to the V data and 



V1555 D31A 



P=0.9167 i=76 




Fig. 3. — y, / lightcurves of eclipsing binaries found in the field M31A. The thin continuous line 
represents for each system the best fit curve (fitted to the V data). 
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V6450 D31A P = 3.076 i = 67 
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Fig. 4. — Continued from Fig.^. 
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Fig. 5. — Continued from Fig.^. 



then for the I data only the zero-point offset was ahowed. 



6.3. Other periodic variables 



For some of the variables preliminary classified as Cepheids (Section 5.3.1), we decided upon 
closer examination to classify them as "other periodic variables" . In Table ^ we present the 
parameters of 6 possible periodic variables other than Cepheids and eclipsing binaries in the 
M31A field, sorted by the increasing period P. For each variable we present its name, 2000.0 
coordinates, value of the variability index J5, period P, error-weighted mean magnitudes V 
and (when available) /. To quantify the amplitude of the variability, we also give the standard 
deviations of the measurements in the V and / bands, cry and cr/. 
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Fig. 6. — V, I lightcurves of Cepheid variables found in the field M31A. The thin continuous line 
represents for each star the best fit Cepheid template (fitted to the V data). 
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21.5 



S2.5 



. 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 . 


- 1 , , , , 1 , , , , 1 , , , , 1 , , , , 1 - 


1 .... 1 .... 1 .... 1 .... 1 ' 



20 

20.5 

21 



0.5 



1.5 3 0.5 

V7466 D31A P = 5.007 A=0.29 



1.5 



- 20.5 




1.5 2 0.5 

V6800 D31A P = 5.077 A=0.35 



22.5 




1.5 2 

V5348 D31A P = 5.346 A=0.41 



21.5 



21.5 




20.5 



20.5 



Fig. 7. — Continued from Fig.^. 



- 14 - 



V8S32 D31A P=5.681 A=0.27 



21 



ai.5 



1 


1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 - 


1 < 1 1 1 1 1 > 1 < 1 > 1 1 1 1 1 < 1 > 1 . 

■ ■ ■ ■ ■ 


■ 1 . 


. . 1 , . . . 1 . 




- 1 , , , , 1 , , , , 1 , , , , 1 , , , , 1 - 



19.5 



0.5 



1.6 3 0.5 

V589 D31A P = 6.192 A = 0.4 



1.5 




21 
21.5 

22 
22.5 

21 



21 



21.5 



22.5 



20.5 






0.5 


1 


1.5 3 
V3142 D31A 




P = 6.244 


0.5 1 
A=0.38 


1.5 2 


. 1 1 1 1 


1 1 1 1 


1 1 1 1 1 1 


1 1 1 1 1 1 1 . 


_ 1 ' ' ' 


1,11.1, 


' ' ' ' 1 ' ' ' ' 1 J 


■/ 




. ? 




my 






■ "" 




; " 






■ 

■ ■ 


■ ■ 

■ 

■ ■ 


r 1 , , , 


. 1 , . 


. . 1 . . . 


, 1 , , , , 1 " 


- , . ■ 

1 . . . 


, 1 .... 1 


■ 

.... 1 .... 1 : 





0.5 


1 


1.5 3 
V2770 D31A 




P = 6.413 


0.5 1 
A=0.29 


1.5 2 


,1.1 


1 , . 1 


1 1 1 1 1 1 


.,11.1, 


1 1 1 1 
1 


.,11.1, 
■ ^ 


1 . 1 1 , . 1 1 1 1 


■ 1 , , , 


. 1 , . 


■ ■ ■ 

. . 1 . . . 


■ ■ 

■■■■T^pii.. 

, 1 .... 1 


■ ■ 

■ ■ 

■ ■ 

1 . . . 


, 1 .... 1 


■ 1. . 

■ ■ 

.... 1 .... 1 . 





0.5 


1 


1.5 2 
V6B42 D31A 




P = 6.482 


0.5 1 
A=0.3 


1.5 2 


1 1 , 1 


1 1 . 1 


1 1 1 1 1 1 


' 1 ' ' ' ' 1 _ 
f 

^if ■ 


1 ' ' ' 


.,11.1, 
■ ■ 


1 . 1 1 , . 1 1 1 1 
■ ■ 










■ j/^ i 

■ 




■ 


1 ■■ 

■ 1 , , , 


. 1 , . 


1 11 

. . 1 . . . 


, 1 . , , . 1 


■ 

1 . . . 


■ ■ 

■ ■ . 

, 1 .... 1 


■ ■ ■ . 

! ■ : 

.... 1 .... 1 ' 





0.5 


1 


1.5 2 
V9473 D31A 




P = 6.582 


0.5 1 
A=0.32 


1.5 2 


- 1 1 , 1 


1 1 . 1 
1 

■ 


1 1 1 1 1 1 


.,11.1,- 
1 


-I'll 


.,11.1, 


1 . 1 1 , . 1 1 1 1 _ 












■ 


■ 


" ■ 

• ■ 

■ 1 , , , 


1 

. 1 , . 


■ ■ 

. . 1 . . . 


, 1 .... 1 


■ ■ 

1 . . . 


, 1 .... 1 


■ ■ 

.... 1 .... 1 ' 





0.5 


1 


1.5 2 





0.5 1 


1.5 2 



20.5 
21 
21.5 
22 



20.5 



21 



21.5 



20 



20.5 



Phase 



Fig. 8. — Continued from Fig.^. 
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Fig. 10. — Continued from Fig.^. 
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Fig. 11. — Continued from Fig.^. 
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Table 1. DIRECT Eclipsing Binaries in M31A 



Name 
(D31A) 


02000.0 

[deg] 


(52000.0 

[deg] 


Js 


P 

[days] 


^max 


Imax 


Ri 


R2 


i 

[deg] 


e 


Comments 


V1555 


11.2717 


41.6462 


1.07 


0.917 


20.64 




0.59 


0.41 


76 


0.00 


V6913 D31B 


V9936 


11.4523 


41.7306 


0.77 


0.930 


21.33 




0.65 


0.34 


61 


0.06 




V4741 


11.3337 


41.7511 


1.43 


1.604 


20.35 


20.16 


0.51 


0.40 


73 


0.03 




V8420 


11.4072 


41.7188 


1.21 


1.931 


20.71 


20.58 


0.42 


0.37 


74 


0.04 




V6024 


11.3611 


41.6841 


0.83 


2.083 


20.93 




0.38 


0.37 


69 


0.00 




V7393 


11.3851 


41.7951 


0.85 


2.786 


20.14 




0.35 


0.31 


80 


0.03 




V6450 


11.3681 


41.7400 


1.39 


3.076 


20.60 


20.63 


0.56 


0.44 


67 


0.04 




V6527 


11.3671 


41.8256 


0.84 


4.180 


20.77 




0.31 


0.31 


76 


0.01 




V5912 


11.3563 


41.7510 


0.88 


5.006 


19.92 




0.67 


0.33 


61 


0.01 




V9840 


11.4495 


41.7034 


0.78 


5.215 


21.51 




0.63 


0.36 


75 


0.03 




V8153 


11.4026 


41.6701 


0.87 


5.794 


21.49 


20.36 


0.56 


0.44 


83 


0.14 




V5407 


11.3464 


41.7504 


1.02 


7.061 


20.33 


19.51 


0.51 


0.44 


57 


0.02 




V538 


11.2420 


41.6695 


0.85 


7.171 


21.07 


20.44 


0.59 


0.40 


66 


0.00 




V4636 


11.3321 


41.7514 


0.80 


8.181 


19.44 


19.43 


0.14 


0.14 


83 


0.18 


DEB 


V6423 


11.3671 


41.7533 


0.85 


11.782 


20.80 


19.89 


0.54 


0.46 


57 


0.15 





Note. — Variable V4636 D31A with period P = 8.181 days is a detached eclipsing binary (DEB) 
with significant ellipticity. V1555 D31A was found as V6913 D31B in Paper I, with identical period, 
Vmax = 20.63 and Imax = 20.06. 
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Table 2. DIRECT Cepheids in M31A 



TVT 


"2000.0 


c 

(>2000.0 




F 




{I) 


A 


Comments 


[u61A) 


[degj 


r 1 ,1 

[degJ 




\days\ 










V3885 


11.3176 


41.7868 


1.13 


3.801 


21.87 


21.63 


0.34 




V4585 


11.3337 


41 6685 


1.18 


3.953 


22.06 


21.71 


0.44 




V8435 


11.4083 


41 6933 


1.05 


4.097 


22.05 


21.37 


0.43 




V3071 


11 3036 


41 6866 


1.05 


4.557 


21.49 




0.23 




V8798 


11.4135 


41.8081 


1.36 


4.861 


22.09 




0.38 




V9833 


11.4461 


41.7914 


1.60 


4.867 


21.22 


20.20 


0.29 




V9531 


11.4361 


41.7111 


1.69 


4.869 


21.55 


20.55 


0.36 




V7466 


11.3867 


41.7830 


1.15 


5.007 


21.42 


20.71 


0.29 




V3584 


11.3125 


41.7674 


0.80 


5.020 


21.09 


20.12 


0.14 




V6800 


11.3780 


41.6476 


0.91 


5.077 


21.82 




0.35 




V5348 


11.3479 


41.6677 


1.65 


5.346 


21.61 


20.75 


0.41 




V8573 


11.4116 


41.7128 


1.35 


5.476 


21.39 


20.50 


0.32 


Ma97 125 


V8232 


11.4011 


41.7727 


1.37 


5.681 


21.17 


19.88 


0.27 


Ma97 124 


V589 


11.2430 


41.6840 


2.42 


6.192 


21.25 


20.25 


0.40 




V3142 


11.3047 


41.7006 


1.83 


6.244 


21.52 


21.09 


0.38 




V2770 


11 2952 

J Lb ^ ty ^ 


41.7592 


1.62 


6.413 


21.25 


20.46 


0.29 




V6842 


11.3770 


41.7058 


1.20 


6.482 


21.65 


20.79 


0.30 




V9473 


11.4337 


41.7366 


1.76 


6.582 


21.16 


20.22 


0.32 


Ma97 127 


V5188 


11.3438 


41.7022 


1.52 


6.776 


21.31 


20.38 


0.23 


Ma97 111 


V1416 


11.2661 


41 6983 


2.98 


6.925 


20.64 


19.57 


0.28 




V3568 


11.3117 


41.7760 


2.10 


7.165 


21.31 


20.73 


0.35 




V5968 


11 3586 


41.7275 


1.64 


8.519 


21.12 


20.05 


0.24 


Ma97 114 


V2242 


11.2854 


41.7508 


1.12 


8.680 


21.47 


20.51 


0.26 




V7523 


11.3917 


41.6581 


1.82 


8.709 


21.05 


20.43 


0.24 


Ma97 121 


V2276 


11 2886 


41 6656 


1.42 


9.803 


20.86 


19.70 


0.17 


V7553 D31B 


7Q1 


Ll.Zl Do 






lU.Uli 




1 Q c;s 
ly.oo 


U.ZD 




V6363 


11.3683 


41.6595 


1.18 


10.593 


21.41 


20.01 


0.26 


Ma97 117 


V4733 


11.3324 


41.7888 


1.99 


10.971 


21.39 


20.11 


0.35 




V9029 


11.4211 


41.7472 


2.95 


11.668 


20.84 


19.71 


0.32 


Ma97 126 


V107 


11.2255 


41.7051 


4.41 


12.525 


20.59 


19.66 


0.37 




V4104 


11.3261 


41.6504 


1.45 


12.557 


21.21 


19.98 


0.23 




V3407 


11.3085 


41.7645 


4.39 


12.801 


20.27 


19.40 


0.39 




V8530 


11.4076 


41.8088 


0.79 


12.809 


21.38 


20.17 


0.22 




V8882 


11.4152 


41.8081 


1.51 


13.097 


21.21 


19.76 


0.35 




V4407 


11.3308 


41.6660 


1.70 


14.112 


19.97 


19.14 


0.14 




V6759 


11.3772 


41.6576 


6.45 


15.479 


20.46 


19.41 


0.51 


Ma97 118 
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6.4. Miscellaneous variables 

In Table Q we present the parameters of 11 miscellaneous variables in the M31A field, sorted 
by the decreasing value of the mean magnitude V. For each variable we present its name, 2000.0 
coordinates, value of the variability index Jsi> 1-2), mean magnitudes V and I. To quantify the 
amplitude of the variability, we also give the standard deviations of the measurements in V and 
/ bands, ay and aj. In the "Comments" column we give a rather broad sub-classification of the 
variability: LP - possible long-period variable (P > 56 days); IRR - irregular variable. Unlike 
the M31B field (Paper I), all of the miscellaneous variables seem to represent the LP type of 
variability. 



6.5. Comparison with other catalogs 

The area of M31A field has not been observed frequently before and the only overlapping 
variable star catalog is given by Magnier et al. (1997, hereafter Ma97). Of 17 variable stars in 
Ma97 which are in our M31A field, we cross-identified 16. Of these 16 stars, one (Ma97 120) we 
did not classify as variable (J^ = 0.22) and one (Ma97 123) was initially classified as variable 
{Js = 1.1), but then failed to classify as a Cepheid. The remaining 14 stars we classified as 
Cepheids (see Table § for cross-identifications) . 

There was also by design a slight overlap between the M31A and M31B fields (Fig.[l8|). Out of 
two eclipsing binaries found in the overlap from the M31A field, V1555 D31A was cross-identified 
as V6913 D31B, with very similar properties of the light curves (see Table ||). V538 D31A turned 
out to have Js = 0.749 in the M31B database, i.e. it very narrowly escaped classification as a 
candidate variable star in this field. There was only one Cepheid from the M31A field in the 
overlap region, V2276 D31A (also Ma97 108), and it was cross-identified as V7553 D31B, again 
with very similar properties of the light curves (see Table |2|). There was a second Cepheid in 
the overlap found in Paper I, V7845 D31B {Js = 0.88), which failed to qualify as a Cepheid 
in the M31A field because of the more stringent requirements for the reduced x^/^dof of the 
fit (Section 5). We also cross- identified three miscellaneous variables (see Table ^), out of four 
detected in the M31A field and five detected in the M31B field, which fell into the overlap region. 



Discussion 



In Fig. 17 we show V,V — I color-magnitude diagrams for the variable stars found in the field 
M31A. The eclipsing binaries and Cepheids are plotted in the left panel and the other periodic 
variables and miscellaneous variables are plotted in the right panel. As expected, some of the 
eclipsing binaries occupy the blue upper main sequence of M31 stars, but there is a group of 
eclipsing binaries with F — / ~ 1.0. The Cepheid variables group near V — I ^ 1.0, with the 
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Table 2 — Continued 



Name 


q;20oo.o 


(^2000.0 


Js 


P 


{V) 


{I) 


A 


Comments 


(D31A) 


[deg] 


[deg] 




[days] 










V5760 


11.3544 


41.7088 


4.49 


16.608 


20.83 


19.69 


0.44 




V5614 


11.3509 


41.7351 


5.43 


20.18 


20.46 


19.37 


0.38 


Ma97 113 


V4452 


11.3319 


41.6506 


1.37 


26.59 


21.56 


19.52 


0.38 




V6165 


11.3632 


41.6997 


5.58 


28.76 


20.05 


18.84 


0.43 


Ma97 116 


V4711 


11.3361 


41.6596 


1.15 


32.29 


21.54 


20.06 


0.38 




V5415 


11.3483 


41.6933 


4.46 


37.06 


20.64 


19.02 


0.48 


Ma97 112 


V9679 


11.4405 


41.7749 


2.46 


42.55 


20.18 


18.75 


0.33 


Ma97 129 



Note. — Variable V2276 D31A (Ma97 108) was found as V7553 D31B in Paper I, with 
period P = 9.482 days, {V) = 20.93 and (/) = 19.77. 



Table 3. DIRECT Other Periodic Variables in M31A 



Name 


02000.0 


(^2000.0 


Js 


P 


V 


/ 


av 




Comments 


(D31A) 


[deg] 


[deg] 




[days] 












V1494 


11.2659 


41.7755 


1.30 


23.15 


21.56 




0.33 






V410 


11.2355 


41.7199 


1.10 


28.62 


20.92 




0.13 






V1911 


11.2800 


41.6786 


1.68 


36.74 


19.99 


19.25 


0.12 


0.14 


RV Tau 


V3368 


11.3104 


41.6779 


2.62 


46.18 


21.14 


20.34 


0.36 


0.23 


RV Tau 


V2977 


11.3018 


41.6835 


1.27 


53.3 


21.21 


20.13 


0.18 


0.28 


RV Tau 


V9659 


11.4421 


41.6874 


1.75 


56.6 


20.77 


19.96 


0.40 


0.27 


RV Tau 
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Table 4. DIRECT Miscellaneous Variables in M31A 



Name 


"2000.0 


^2000.0 


Js 


V 








Comments 


(D31A) 


[deg] 


[deg] 














V3222 


11.3081 


41.6494 


1.40 


19.11 


17.21 


0.09 


0.04 


V8123 D31B 


V3901 


11.3216 


41.6633 


1.61 


19.25 


16.72 


0.10 


0.05 


LP 


V2109 


11.2844 


41.6975 


1.67 


20.07 


19.13 


0.11 


0.07 


LP 


V7718 


11.3926 


41.7313 


2.02 


20.23 


19.59 


0.15 


0.10 


LP 


V8415 


11.4050 


41.7867 


1.98 


20.37 


19.65 


0.17 


0.12 


LP 


V5764 


11.3563 


41.6507 


1.36 


20.75 


20.20 


0.14 


0.12 


LP 


V2570 


11.2912 


41.7629 


1.82 


21.42 


19.60 


0.53 


0.20 


LP 


V541 


11.2424 


41.6594 


2.47 


21.51 


19.20 


0.54 


0.13 


V5897 D31B 


Vll 


11.2240 


41.6504 


1.27 


21.72 


19.31 


0.47 


0.12 


V5075 D31B 


V9459 


11.4322 


41.7689 


1.71 


21.84 


19.23 


0.54 


0.18 


LP 


V1032 


11.2565 


41.6655 


1.72 


21.88 


19.43 


0.43 


0.14 


LP 



Note. — 



Variables V3222, V541 and Vll were also found in Paper L 
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Fig. 14. — V, I lightcurves of other periodic variables found in the field M31A. 
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Fig. 15. — V, I lightcurves of miscellaneous variables found in the field M31A. 
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Fig. 17. — V, V — I color-magnitude diagrams for the variable stars found in the field M31A. 
The eclipsing binaries and Cepheids are plotted in the left panel and the other periodic variables 
and miscellaneous variables are plotted in the right panel. The dashed line corresponds to the I 
detection limit of / ~ 21 mag. 



exception of possibly highly reddened system V4452 D31A. The other periodic variable stars 
have positions on the CMD similar to the Cepheids. The miscellaneous variables are scattered 
throughout the CMD and might represent few classes of variability, but most of them are red with 
V — I = 1.8 — 2.6, and are probably Mira variables. 

In Fig.|l^ we plot the location of eclipsing binaries and Cepheids in the fields M31A and 
M31B, along with the blue stars {B — V < 0.4) selected from the photometric survey of M31 by 
Magnier et al. (1992) and Haiman et al. (1994). The sizes of the circles representing the Cepheids 
variables are proportional to the logarithm of their period. As could have been expected, both 
types of variables group along the spiral arms, as they represent relatively young populations of 
stars. However, in the field M31A there is a group of Cepheids located outside the population of 
the blue stars (at ~ 11.3,41.76 deg and below). Also, these Cepheids as well as the Cepheids 
located in the inner spiral arm of M31 (field M31B at RA~ 11.1 deg) have on average shorter 
periods than the Cepheids located in the outer spiral arm. We will explore various properties of 
our sample of Cepheids in the future paper (Sasselov et al. 1998, in preparation). 



We would like to thank the TAC of the Michigan-Dartmouth-MIT (MDM) Observatory 
and the TAC of the F. L. Whipple Observatory (FLWO) for the generous amounts of telescope 
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Fig. 18. — Location of eclipsing binaries (filled squares) and Cepheids (open circles) in the fields 
M31A and M31B, along with the blue stars {B — V < 0.4) selected from the photometric survey of 
M31 by Magnier et al. (1992) and Haiman et al. (1994). The sizes of the circles representing the 
Cepheids variables are proportional to the logarithm of their period. 
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